Modeling  and  Reduction  of  Shocks  on  Electronic 
Components  Within  a  Projectile 


by  Vinod  Chakka,  Mohamed  B.  Trabia,  Brendan  O’Toole, 
Srujanbabu  Sridharala,  Samaan  Ladkany,  and  Mostaflz  Chowdhury 


ARL-RP-217 


August  2008 


A  reprint  from  the  International  Journal  of  Impact  Engineering,  vol.  35,  pp.  1326-1338,  2008. 


Approved  for  public  release;  distribution  is  unlimited. 


NOTICES 

Disclaimers 

The  findings  in  this  report  are  not  to  be  construed  as  an  official  Department  of  the  Army  position  unless 
so  designated  by  other  authorized  documents. 

Citation  of  manufacturer’s  or  trade  names  does  not  constitute  an  official  endorsement  or  approval  of  the 
use  thereof 


Destroy  this  report  when  it  is  no  longer  needed.  Do  not  return  it  to  the  originator. 


Army  Research  Laboratory 

Adelphi,MD  20783-1197 


ARL-RP-217 


August  2008 


Modeling  and  Reduction  of  Shocks  on  Electronic 
Components  Within  a  Projectile 

Vinod  Chakka,  Mohamed  B.  Trabia,  Brendan  O’Toole, 
Srujanbabu  Sridharala,  and  Samaan  Ladkany 
University  of  Nevada 

Mostafiz  Chowdhury 

Weapons  and  Materials  Research  Directorate,  ARL 


A  reprint  from  the  International  Journal  of  Impact  Engineering,  vol.  35,  pp.  1326-1338,  2008. 


Approved  for  public  release;  distribution  is  unlimited. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering 
and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  the  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson 
Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to 
comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY) 

August  2008 


2.  REPORT  TYPE 

Reprint 


4.  TITLE  AND  SUBTITLE 

Modeling  and  Reduction  of  Shocks  on  Electronic  Components  Within  a  Projectile 


3.  DATES  COVERED  (From  -  To) 

January  2005-January  2007 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 

Vinod  Chakka,  Mohamed  B.  Trabia,  Brendan  O’Toole,  Srujanbabu  Sridharala, 
Samaan  Ladkany,^  and  Mostafiz  Chowdhury 


5d.  PROJECT  NUMBER 

DAAD19-03-2-0007 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Research  Laboratory 
ATTN:  AMSRD-ARL-WM-MB 
Adelphi,MD  20783-1197 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

ARL-RP-217 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 


11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 


13.  SUPPLEMENTARY  NOTES 

A  reprint  from  the  International  Journal  of  Impact  Engineering,  vol.  35,  pp.  1326-1338,  2008. 
*Department  of  Mechanical  Engineering,  University  of  Nevada,  Las  Vegas,  NV  89154-4027 
^Department  of  Civil  and  Environmental  Engineering,  University  of  Nevada,  Las  Vegas,  NV  89154-4027 


14.  ABSTRACT 

Electronic  components  within  a  projectile  are  subjected  to  severe  loads  over  an  extremely  short  duration  during  the  launch 
process.  Failure  of  these  components  during  launch  can  result  in  negative  effects  on  the  mission  of  the  projectile.  While 
experimental  data  can  be  helpful  in  understanding  failure  of  electronic  components  within  a  projectile,  collecting  such  data  are 
[sic]  usually  difficult.  There  are  also  limitations  on  the  reliability  of  sensors  under  these  circumstances.  Finite  element 
modeling  (FEM)  can  offer  a  means  to  better  understand  the  behavior  of  these  components.  It  can  also  be  used  to  develop  better 
shock  mitigation  features  into  the  projectile  design.  This  research  has  two  objectives.  The  first  objective  is  to  develop  an  FEM 
that  one  describes  the  interaction  of  a  typical  projectile  with  the  gun  barrel  during  launch.  The  projectile  includes  a  payload  of  a 
one-pound  mass  representing  a  typical  electronic  package  supported  by  a  plate.  The  second  objective  of  this  work  is  to 
investigate  the  use  of  composite  plates  to  support  electronic  payload  as  a  means  to  reduce  the  transmitted  shocks  during  the 
projectile  launch  event.  The  proposed  plate  has  carbon  fibers  embedded  in  an  epoxy  matrix.  A  parametric  study  of  the  effects 
of  varying  the  thickness  of  the  supporting  plate  and  the  fiber  volume  fraction  on  the  accelerations  and  stresses  is  included. 

Results  of  the  study  are  used  to  reach  general  recommendations  regarding  reducing  failure  of  electronic  components  within  a 
projectile. _ 


15.  SUBJECT  TERMS 

projectile,  finite  element  modeling,  electronic  component,  shock,  composites 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT 

18.  NUMBER 

OF  PAGES 

19a.  NAME  OF  RESPONSIBLE  PERSON 

Mostafiz  R.  Chowdhury 

a.  REPORT 

UNCLASSIFIED 

b.  ABSTRACT 

UNCLASSIFIED 

c.  THIS  PAGE 

UNCLASSIFIED 

UL 

24 

19b.  TELEPHONE  NUMBER  {Include  area  code) 

(301)  394-6308 

standard  Form  298  (Rev.  8/98) 
Prescribed  by  ANSI  Std.  Z39.18 


Available  online  at  www.sciencedirect.com 


ScienceDirect 


International  Journal  of  Impact  Engineering  35  (2008)  1326-1338 


INTERNATIONAL 
JOURNAL  OF 


IMPACT 

ENGINEERING 


www.elsevier.com/locate/ijimpeng 


Modeling  and  reduction  of  shocks  on  electronic  components 

within  a  projectile 

Vinod  Chakka^,  Mohamed  B.  Trabia^’*,  Brendan  O’Toole^,  Srujanbabu  Sridharala^, 

Samaan  Ladkany’’,  Mostafiz  Chowdhury^^ 

^Department  of  Mechanieal  Engineering,  University  of  Nevada,  Las  Vegas,  Las  Vegas,  NV  89154-4027,  USA 
^Department  of  Civil  and  Environmental  Engineering,  University  of  Nevada,  Las  Vegas,  Las  Vegas,  NV  89154-4027,  USA 
^AMSRL-WM-MB  (ALC),  Ordnance  Materials  Branch,  Weapons  &  Materials  Research  Directorate,  US  Army  Research  Laboratory, 

2800  Powder  Mill  Road,  Adelphi,  MD  20783-1145,  USA 

Received  16  October  2006;  received  in  revised  form  12  January  2007;  accepted  24  July  2007 
Available  online  3  August  2007 


Abstract 

Electronic  components  within  a  projectile  are  subjected  to  severe  loads  over  an  extremely  short  duration  during  the  launch  process. 
Failure  of  these  components  during  launch  can  result  in  negative  effects  on  the  mission  of  the  projectile.  While  experimental  data  can  be 
helpful  in  understanding  failure  of  electronic  components  within  a  projectile,  collecting  such  data  are  usually  difficult.  There  are  also 
limitations  on  the  reliability  of  sensors  under  these  circumstances.  Finite  element  modeling  (FEM)  can  offer  a  means  to  better  understand 
the  behavior  of  these  components.  It  can  also  be  used  to  develop  better  shock  mitigation  features  into  the  projectile  design.  This  research 
has  two  objectives.  The  first  objective  is  to  develop  an  FEM  that  one  describes  the  interaction  of  a  typical  projectile  with  the  gun  barrel 
during  launch.  The  projectile  includes  a  payload  of  a  one-pound  mass  representing  a  typical  electronic  package  supported  by  a  plate.  The 
second  objective  of  this  work  is  to  investigate  the  use  of  composite  plates  to  support  electronic  payload  as  a  means  to  reduce  the 
transmitted  shocks  during  the  projectile  launch  event.  The  proposed  plate  has  carbon  fibers  embedded  in  an  epoxy  matrix.  A  parametric 
study  of  the  effects  of  varying  the  thickness  of  the  supporting  plate  and  the  fiber  volume  fraction  on  the  accelerations  and  stresses  is 
included.  Results  of  the  study  are  used  to  reach  general  recommendations  regarding  reducing  failure  of  electronic  components  within  a 
projectile. 

(C)  2007  Elsevier  Etd.  All  rights  reserved. 

Keywords:  Projectile;  Modeling;  Shock;  Finite  element;  Electronic  components 


1.  Introduction 

During  the  last  thirty  years,  the  US  Army  has  been 
developing  “smart  artillery”  munitions,  which  contain 
sophisticated  embedded  electronic  systems.  Unfortunately, 
the  artillery  environment  is  extremely  harsh  as  the 
munitions  must  operate  in  temperatures  ranging  from 
— 60°F  to  160°F.  The  projectiles  are  subjected  to  quasi¬ 
static  axial  loads  in  excess  of  15,000  g  during  the  launch 
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process  which  is  augmented  by  transient  loads  of  up  to 
5000  g  at  the  muzzle  exit  [1].  Projectiles  can  also  spin  at  up 
to  300  revolutions  per  second.  As  a  projectile  travels  down 
the  gun  barrel,  it  may  also  be  subject  to  off-axis  loads  from 
impacts  with  the  gun  tube  walls  caused  by  balloting  [2]. 

Heaslip  and  Punch  [3]  concluded  that  there  is  consider¬ 
able  evidence  to  suggest  that  a  large  percentage  of  portable 
electronic  products  fail  due  to  impact  or  shock  during  use. 
Some  of  the  electronic  component  failure  problems  can  be 
reduced  by  enclosing  the  components  within  foam.  This 
approach,  however,  reduces  the  requirement  that  the 
electronic  components  should  also  be  upgradeable  and 
replaceable  without  swapping  an  entire  subsystem  of  the 
projectile.  These  challenges  present  significant  problems  for 
designers  who  typically  resort  to  the  use  of  numerical 
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simulations  to  provide  guidance  on  these  complex  design 
issues. 

One  of  the  difficulties  facing  the  US  Army  and  its 
contractors  is  identifying  the  correct  gun  launch  loads  to 
electronic  component  manufacturers  prior  to  the  final 
design  of  the  projectile.  In  the  past,  pressure-time  curves 
and  peak  acceleration  values  were  provided  to  contractors. 
Contractors  used  the  peak  values  to  perform  static  analysis 
and  quasi-static  in  centrifuge  tests.  The  dynamics  of  the 
projectile  structure  [3],  particularly  during  the  muzzle  exit 
transient,  were  usually  neglected.  As  a  result,  programs  like 
the  US  Army’s  Excalibur  and  SAD  ARM  experienced 
numerous  failures  of  sensitive  equipment  during  the  early 
stages  of  development.  Failure  of  several  sensitive  compo¬ 
nents  was  traced  to  the  muzzle  exit  event  in  the  Excalibur 
program  [4]. 

Few  researchers  considered  creating  a  fully  detailed 
model  of  a  projectile  launch  event.  This  research  has  two 
objectives: 

-  Develop  a  gun-launched  projectile  finite  element  model 
(FEM)  capable  of  predicting  shocks  transmitted  to  on¬ 
board  electronic  components  during  the  launch  process. 
The  projectile  is  subjected  to  a  realistic  launch  pressur¬ 
e-time  history.  The  model  includes  friction  between  the 
gun  barrel  inner  surface  and  the  projectile  as  well  as  the 
flexibility  of  the  gun  barrel.  The  effects  of  these  factors 
on  the  vibrations  of  the  electronic  package  are  also 
considered. 

-  Investigate  the  effects  of  composite  laminate  parameters 
on  reducing  the  transmitted  shocks  to  embedded 
electronic  components  through  passive  damping. 

The  following  is  a  review  of  relevant  research.  Simkins  et 
al.  [5]  studied  dynamics  of  gun  tubes,  especially  interaction 
between  strains  in  the  wall  of  the  gun  tube  and  projectile 
velocity  to  predict  resonant  conditions.  Hopkins  and 
Wilkerson  [6]  examined  the  results  of  a  series  of  experi¬ 
ments  conducted  in  an  attempt  to  reduce  the  dynamic 
motions  of  the  M256  gun  system  during  firing.  Data 
collected  during  these  experiments  included  the  motion  of 
the  gun  tube  and  breech  mechanism.  An  axisymmetric 
transient  FEM  of  the  entire  system  was  developed  and  used 
to  explain  the  experimental  results.  Hollis  [7]  developed  a 
two-dimensional  quasi-static  model  of  a  training  projectile. 
The  projectile  was  redesigned  to  reduce  stresses.  Cordes  et 
al.  [2]  presented  a  simplified  model  of  a  projectile  using 
shell  elements.  Natural  frequency  effects  were  considered  in 
this  model.  Kessler  and  Spearing  [8]  studied  an  autono¬ 
mous  flying  vehicle  that  is  launched  from  an  artillery  shell, 
and  then  deployed  over  a  battlefield  to  capture  images. 
They  developed  the  structural  design  and  testing  of  the  aft 
section  of  this  vehicle,  which  is  subjected  to  high  impulsive 
inertial  loads.  A  quasi-dynamic  FEM  of  the  aft  section  was 
developed  whose  results  were  compared  with  launch  data. 
Chowdhury  et  al.  [1]  successfully  used  a  3D  FE  gun-launch 
simulation  model — including  balloting,  spinning,  and  gun 


mount  positioning — of  a  surrogate  Excalibur  projectile  in 
support  of  its  design  development.  In  their  paper,  the 
impacts  of  an  FEM  in  locating  trouble-free  mounting 
locations  for  sensitive  electronic  components  and  the  need 
for  parametric  study  in  identifying  sensitive  factors 
affecting  the  muzzle  exit  dynamics  of  projectile  was 
demonstrated.  Petersen  [9]  proposed  a  test  pan  to  gather 
information  on  the  dynamics  of  the  155-mm  Advanced 
Gun  System  barrel  as  a  projectile  travels  down  the  gun  tube 
to  corroborate  modeling  efforts  using  sensors  on  the  gun 
barrel  to  measure  the  position,  acceleration,  and  strain 
during  firing.  The  acquired  data  are  compared  to  barrel- 
projectile  interaction  FEM.  Recently,  an  analytic  approach 
is  proposed  to  investigate  the  dynamic  behavior  of 
laboratory  rail  guns  resulting  from  launching  a  projectile 
[10].  The  rail  gun  is  modeled  as  a  beam  with  an  elastic 
foundation  with  cantilevered  support  at  the  breech  end  of 
rails.  The  analysis  showed  that  structural  response  of  the 
rail  is  governed  by  a  transient  fourth-order  differential 
equation.  Flexibility  of  the  gun  was  not  considered  in  most 
of  these  models. 

In  addition  to  developing  an  accurate  model  of  the 
projectile  and  gun  barrel,  it  is  important  to  explore  means 
of  reducing  transmitted  vibrations  to  electronic  payloads. 
Most  of  the  research  in  this  area  concerns  the  use  of  smart 
materials  to  passively  dampen  vibrations  such  as  dynamic 
vibration  absorbers  or  resonant  circuit  shunted  piezo¬ 
ceramics  [11].  On  the  other  hand,  few  researchers 
considered  modifying  characteristics  of  composite  struc¬ 
ture  to  passively  reduce  vibrations.  For  example,  Trindade 
[12]  used  genetic  algorithms  to  optimize  the  geometry  of  a 
laminated  composite  beam  for  passive  damping.  None  of 
the  surveyed  literature  addressed  the  problem  of  shock 
loading  within  a  projectile  launch  event. 

The  remainder  of  this  paper  is  divided  as  follows:  Section 
2  presents  description  of  the  projectile  and  gun.  An  FEM 
for  the  projectile  and  gun  is  developed  in  Section  3.  An 
approach  for  the  design  of  the  electronic  payload 
supporting  composite  plate  is  presented  in  Section  4. 
Conclusions  and  recommendations  for  future  work  are 
listed  in  Section  5. 

2.  Description  of  the  projectile  and  gun 

The  US  Army  ARDEC  Tank-automotive  &  Armaments 
Command-Armaments  Research,  Development  &  Engi¬ 
neering  Center  (TACOM-ARDEC)  provided  the  data  for 
the  projectile  under  consideration,  using  an  instrumented 
projectile.  Fig.  1(a),  to  collect  data  within  and  outside  the 
gun  including  pressure  and  acceleration  time  histories.  The 
projectile  is  composed  of  several  components:  M795  Body, 
ogive,  nacelle,  and  windshield.  All  these  components  are 
threaded  into  each  other  or  bolted  together.  A  band  on  the 
M795  maintains  contact  between  the  projectile  and  the 
walls  of  the  gun.  The  properties  for  all  parts  are  listed  in 
Tables  1  and  2.  The  actual  projectile  has  various  electronic 
components  occupying  its  top  half  in  addition  to  pressure 
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Fig.  1.  Sectional  view  of  the  original  projectile  and  its  simplified  model. 


M795  Body 


Table  1 

Material  properties  of  projectile  parts  and  gun  barrel 


Part  name 

Material 

Specific  weight 
(Ibf/in^*) 

Wind  shield 

Ultem  2300  (30% 
glass) 

0.055 

Nacelle 

Aluminum 

7075-T6511 

0.103 

Ogive,  M795,  plate, 
payload,  band,  gun 
barrel 

4030  Steel 

0.283 

transducers  at  the  bottom  and  mid-section  of  the  projectile 
to  measure  axial  and  radial  pressures.  The  modeling  was 
simplified  by  eliminating  the  exact  details  of  all  the 
electronic  components  and  focusing  on  the  response  of  a 
single  generic  electronic  component.  Based  on  suggestions 
of  TACOM-ARDEC,  it  was  decided  to  eliminate  these 
electronic  components  and  to  adjust  the  mass  of  nacelle 
accordingly.  A  one-pound  mass,  which  corresponds  to  a 
typical  electronic  package,  is  added  as  shown  in  Fig.  1(b). 
This  payload  is  supported  by  a  steel  plate  that  is  fully 
attached  to  the  Ogive  (Fig.  1(b)). The  thickness  of  this  plate 
is  similar  to  what  is  typically  used  in  industry.  Critical 
dimensions  of  the  projectile  are  shown  in  Fig.  2.  The  inner 


diameter  of  the  Ogive  is  slightly  reduced  to  provide  an 
attachment  point  for  the  plate.  The  maximum  diameter  of 
the  projectile  is  6.22  in.  at  the  band,  which  is  equal  to  the 
inner  diameter  of  the  gun  barrel.  The  total  height  of  the 
projectile  is  20.5  in.  The  gun  barrel  is  modeled  as  a  long 
hollow  cylindrical  tube.  The  inside  and  outside  diameters 
of  the  gun  barrel  are  6.22  and  11. 62  in,  respectively.  The 
total  length  of  the  gun  barrel  is  264  in.  The  gun  barrel  is 
supported  by  two  bearings  at  10  and  40  in.  from  the  bottom 
of  the  gun  barrel  (breech  of  the  gun)  as  shown  in  Fig.  3. 
The  gun  barrel  model  used  in  the  simulation  is  a 
simplification  of  the  actual  gun  tube  considering  that  the 
bore  evacuator  has  an  effect  on  launch  dynamics  [13]. 

3.  Development  of  the  FEM 

3.1.  Meshing 

LS-Dyna  [14]  is  used  to  simulate  the  motion  of  the 
projectile.  Most  components  of  the  projectile,  except  for 
the  base  of  the  M795  Body,  the  payload,  and  the  support 
plate,  are  modeled  by  creating  a  two-dimensional  sketch 
for  each  of  them.  Each  sketch  is  meshed  using  quadrilateral 
elements.  These  sketches  are  then  swept  360°  to  create 
corresponding  volumes,  which  are  meshed  using  eight-node 
brick  elements.  The  two-dimensional  mesh  is  then  deleted. 
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Table  2 

Mechanical  characteristics  of  projectile  parts  and  gun  barrel 


Part  name 

Young’s  modulus  (psi) 

Poisson’s  ratio 

Yield  stress  (psi) 

Tangent  modulus  (psi) 

Wind  shield 

8.00E  +  05 

0.40 

24.5E  +  03 

- 

Nacelle 

1.04E  +  07 

0.33 

68.0E  +  03 

1.85E  +  04 

Ogive,  M795,  plate,  payload,  band,  gun  barrel 

2.90E  +  07 

0.32 

120E  +  03 

5.21E  +  04 

00 

ho 

G) 


o 

GO 

CD 


Fig.  2.  Sectional  view  of  the  meshed  projectile  model. 


Fig.  3.  Meshed  gun  barrel  with  boundary  conditions  applied. 


Nodes  are  merged  at  the  beginning  and  end  of  the  swept 
volume  to  eliminate  duplicate  nodes.  A  sectional  view  of 
the  meshed  projectile  is  shown  in  Fig.  2. 


For  the  gun  tube  model,  a  two-dimensional  mesh  is 
created  at  the  base  of  the  gun.  This  mesh  is  extruded  along 
the  axis  of  the  gun  barrel.  About  124  layers  of  eight-node 
brick  elements  are  created  in  this  way.  A  similar  approach 
is  used  to  create  the  base  of  the  M95  Body,  the  payload, 
and  the  support  plate.  Fig.  3  depicts  a  meshed  gun  barrel. 
All  projectile  components  and  gun  models  have  40 
elements  in  the  circumferential  direction. 

3.2.  Loads  and  boundary  conditions 

The  ballistics  of  the  gun  barrel  is  represented  by  a  high-g 
impulsive  thrust  driving  the  base  of  the  projectile.  Fig.  4 
shows  a  typical  pressure-time  history  that  was  measured  at 
the  base  of  the  projectile  using  the  sensor  shown  in  Fig.  1 
[1].  Base  pressure  distributions  on  the  model  are  shown  in 
Figs.  5  and  6.  Gravity  load  is  applied  in  the  negative  x-axis 
of  the  global  frame.  Fig.  7. 

Interfacing  between  various  contact  surfaces  within  the 
projectile  is  defined  by  selecting  a  set  of  elements  on  each  of 
the  two  contacting  parts  that  are  adjacent  to  the  contact 
area.  The  type  of  contact  and  the  role  acted  by  each 
part  are  then  introduced.  Table  3  summarizes  contact 
definitions  between  various  parts  of  the  model.  The 
coefficient  of  friction  between  the  gun  barrel  and  the  band 
is  set  equal  to  0.1. 

As  described  in  the  previous  section,  the  gun  barrel  is 
supported  by  two  bearings  at  10  and  40  in.  from  the  breech 
of  the  gun  barrel.  These  boundary  conditions  are  simulated 
by  identifying  a  band  of  elements  at  each  of  these  two 
locations.  All  degrees  of  freedom  for  nodes  on  the  outer 
surface  of  these  two  bands.  Fig.  3,  are  constrained. 

3.3.  Analysis  approach 

The  projectile  launch  event  is  modeled  in  two  phases.  In 
the  first  phase,  only  the  gun  barrel  is  allowed  to  deflect  due 
to  gravity  in  the  negative  x-direction.  This  analysis  is 
conducted  for  0.5  s,  which  allows  the  gun  to  reach  its  static 
deflection  position.  Fig.  7  shows  the  original  and  deflected 
gun  barrels. 

In  the  second  phase  of  the  analysis,  the  pressure  curve  of 
Fig.  4  is  applied  in  a  normal  direction  to  all  the  surfaces  of 
the  base  portion  of  the  projectile  below  the  band  as  shown 
in  Figs.  5  and  6.  The  initial  conditions  of  this  phase  are 
obtained  from  the  results  of  the  first  phase  of  the  analysis. 
Gravity  load  continues  to  be  active  in  this  phase  also.  The 
duration  of  the  simulation  for  the  second  phase  is  0.02  s. 
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Fig.  5.  Pressure  loads  on  the  base  of  M795  and  the  band. 


Fig.  6.  Pressure  loads  on  the  slanted  surfaces  of  M795. 


The  pressure  pushes  the  projectile  through  the  deformed 
gun  barrel  for  the  first  0.0125  s  of  this  phase.  The  gun  is  not 
internally  pressurized  as  the  effect  of  the  pressure  is  found 
to  be  minimal.  The  initial  position  of  the  projectile  is 
chosen  to  ensure  that  the  lower  edge  of  the  band  is  aligned 
with  the  tip  of  the  gun  barrel  when  the  pressure  curve  is 
first  equal  to  zero  (time  0.0125  s  after  the  beginning  of  the 
second  phase).  Several  trials  result  in  the  conclusion  that 
the  bottom  of  the  projectile  should  be  placed  at  73  in.  from 
the  breech  of  the  gun  barrel.  Results  of  the  FEM  confirm 
this  initial  location.  Fig.  8  shows  the  deformed  shape  of  the 
gun  barrel  in  the  x-direction  (gravity  direction)  at  various 
time  instants.  The  figure  shows  that  gravity  effects 
dominate  the  deformation.  However,  as  the  projectile 
moves  inside  the  gun  barrel  it  slightly  straightens  the  gun. 

Acceleration,  velocity,  and  displacement  data  are  re¬ 
corded  for  nodes  on  the  payload  and  the  nacelle  as  shown 


in  Fig.  9.  Fig.  10  shows  typical  experimental  (unfiltered) 
acceleration  results  for  projectile  of  the  type  used  in  this 
study.  Integrating  this  curve  yields  velocity  and  displace¬ 
ment  curves  that  are  close  to  those  shown  in  Figs.  11  and 
12,  respectively.  However,  acceleration  results.  Figs.  13  and 
14,  have  significantly  higher  level  of  noise  than  what  is 
observed  in  the  experimental  results.  Increasing  the  number 
of  elements  in  the  axial,  tangential,  and  radial  directions 
was  explored  as  a  means  to  reduce  this  acceleration  noise  in 
the  FEM  results.  Such  approach  however  produced  limited 
improvement  in  the  results,  while  proving  to  be  computa¬ 
tionally  prohibitive. 

Acceleration  results  are  therefore  filtered.  A  modal 
analysis  of  the  projectile  is  conducted.  Modal  analysis 
results  are  listed  in  Table  4.  Fig.  15  shows  the  mode  shape 
of  the  first  bending  frequency.  It  is  decided  to  filter 
acceleration  results  at  a  frequency  of  6000  Hz.  This  value  is 
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Contour  (Analysis  system) 
Displacement  (Mag) 

— 2.748E-01 

- 2.442E-01 

-2.137E-01 
-1  832E-01 
-1  526E-01 
-1  221E-01 
-9158E-02 
-6106E-02 
-3  053E-02 

- O.OOOE+00 

■  No  result 
Max=  2.748E-01 
Min  =  0  000E>00 


X 

Z  ^  Y 


Fig.  7.  Static  deflection  of  the  gun  barrel. 


Table  3 

Contact  characteristics 


Master 

Slave 

Contact  type 

Nacelle 

Windshield 

Tied  surface  to  surface 

Ogive 

Nacelle 

Tied  surface  to  surface 

Plate 

Payload 

Tied  surface  to  surface 

Ogive 

Plate 

Tied  surface  to  surface 

M795  body 

Ogive 

Tied  surface  to  surface 

Gun  barrel 

Band 

Surface  to  surface 

Fig.  8.  Deformed  shape  of  the  gun  barrel  in  the  x-direction. 


Node  on  Projectile 


Node  on  Payload 


chosen  as  it  is  larger  than  all  major  frequencies  of  the 
projectile.  Figs.  16  and  17  show  the  filtered  accelerations 
for  the  same  nodes  of  Figs.  13  and  14,  respectively.  These 
results  are  comparable  to  those  observed  in  test  firing  of 
projectiles.  While  accelerations  of  these  two  nodes  are  of 
the  same  order  of  magnitude,  the  payload  exhibits 
significantly  lower-order  frequencies  due  to  the  difference 


Fig.  9.  Locations  of  data  collection. 


in  stiffness  between  the  projectile  and  the  supporting  plate. 
The  results  of  Fig.  17  also  show  that  the  payload 
experiences  accelerations  of  relatively  higher  magnitude 
after  the  projectile  exits  the  muzzle.  This  “ringing” 


1332 


V.  Chakka  et  al  /  International  Journal  of  Impact  Engineering  35  (2008)  1326-1338 


phenomenon  is  of  special  concern  to  failure  of  electronic 
components  in  this  phase. 

4.  Design  of  the  electronic  payload  supporting  composite 
plate 

The  results  of  the  previous  section  show  that  accelera¬ 
tions  of  the  electronic  payload  can  be  divided  into  two 
phases:  inside  the  gun  barrel  and  after  muzzle  exit.  In  the 
first  phase,  accelerations  are  governed  by  the  input 
pressure  curve  and  cannot  be  seriously  affected  by 


choosing  a  different  supporting  plate.  In  the  second 
phase,  the  projectile  and  the  electronic  payload  experience 
free  vibration  as  the  pressure  load  rapidly  disappears. 
The  associated  oscillatory  accelerations  often  lead  to  the 
failure  of  these  components.  Therefore,  choosing  the 
appropriate  plate  characteristics  affects  the  vibration. 
This  is  especially  important  since  experimental  results 
reported  that  failure  of  electronic  components  usually 
occurs  immediately  after  the  projectile  leaves  the  gun 
barrel.  A  desirable  plate  performance  should  combine 
low  peak  accelerations  and  low  frequency  after  the 
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Fig.  13.  Axial  acceleration  for  a  node  on  the  nacelle/projectile. 


projectile  leaves  the  gun  to  ensure  the  survivability  of  the 
electronic  components  it  carries. 

A  parametric  study  to  explore  the  effect  of  using  a 
composite  material  plate  on  reducing  the  accelerations 
transmitted  to  the  payload  is  presented  in  this  section.  The 
plate  is  made  of  standard  modulus  carbon  fibers  embedded 
in  an  epoxy  matrix.  Three  layers  are  used  to  create  the 
laminate.  Individual  fiber  and  epoxy  properties  used  for  the 
analysis  are  listed  in  Table  5.  The  plate  is  designed  as  a 
quasi-isotropic  composite  laminate  for  which  material 
properties  are  isotropic  in  plane.  It  is  decided  to  use  a 
three-laminae  composite  for  the  supporting  plate.  The  fiber 


orientation  angle  of  the  laminae  of  9°,  0°,  and  —0°,  Fig.  18, 
is  suggested  in  [15]  for  this  type  of  application.  Exploded 
and  stacked  views  of  the  three  laminae  are  shown  in  Fig.  18. 
Appendix  A  lists  the  equations  used  to  calculate  the 
properties  of  the  composite  laminate.  Solid  brick  elements 
are  used  to  model  the  composite  plate.  The  fiber  angle 
orientation  of  60°,  0°,  and  —60°  is  used  in  this  work. 
Ongoing  experimental  work  shows  that  the  characteristics  of 
the  supporting  plate  material  are  not  sensitive  to  strain  rate 
variations  within  the  current  projectile  environment  [16]. 

Initially,  a  composite  plate  having  30%  fiber  volume 
fraction  and  0.2  in.  plate  thickness,  which  is  the  same 
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1.50E+05 


1.00E+05 


5.00E+04 


2  O.OOE+00 
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< 


-5.00E+04 


-1.00E+05 


-1.50E+05 


E-01 


Time  (sec) 

Fig.  14.  Axial  acceleration  for  a  node  on  the  payload. 


Table  4 

Modal  results  of  the  projectile 


Mode 

number 

Frequency 

(Hz) 

Type  of  mode 

1 

1786 

1st  Bending  mode  of  payload  and  plate 

2 

1863 

2nd  Bending  mode  of  payload  and  plate 

3 

2008 

3rd  Bending  mode  of  payload  and  plate 

4 

3145 

1st  Bending  mode  of  the  projectile 

5 

3711 

1st  Torsional  mode  of  projectile 

6 

5211 

1st  Axial  mode  of  the  projectile 

7 

5469 

1st  Radial  mode  of  the  projectile 

CO 

"O) 


2 

8 

O 

< 


20E-01 


O.OOE+00 
5.0I 

-5.00E+03 
-1.00E+04 

Time  (sec) 

Fig.  16.  Filtered  axial  acceleration  for  a  node  on  the  nacelle/projectile. 


2  5.00E+03 
0 

8  O.OOE+00 
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-1.00E+04 


E-01 


Time  (sec) 

Fig.  17.  Filtered  axial  acceleration  for  a  node  on  the  payload. 


Fig.  15.  Mode  shape  of  the  first  bending  frequency  of  the  projectile. 

thickness  of  the  original  steel  plate,  is  used.  The  accelera¬ 
tion  results,  Fig.  19,  show  a  significant  reduction  of 
magnitude  and  frequencies  of  the  axial  accelerations  for 
the  node  on  the  payload.  As  expected,  the  acceleration  of 


the  node  on  the  projectile  body  is  practically  unchanged 
when  compared  to  the  case  of  a  projectile  with  a  steel 
supporting  plate. 

The  effect  of  varying  the  fiber  volume  fraction  and  plate 
thickness  on  the  transmitted  accelerations  is  explored  in  the 
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Table  5 

Properties  of  fiber  and  matrix 


Young’s  modulus  of  fiber  (Ef) 

3.19E  +  07psi 

Young’s  modulus  of  matrix  (E^f) 

5.22E  +  05psi 

Poisson’s  ratio  of  fiber  (vf) 

0.2 

Poisson’s  ratio  of  matrix  (vm) 

0.35 

Shear  modulus  of  fiber  (Gf) 

1.32E  +  07psi 

Shear  modulus  of  matrix  (Gm) 

1.92E  +  05psi 

Density  of  the  fiber  (pf) 

1.65E-41b/in^ 

Density  of  the  matrix  (pm) 

1.12E-41b/m^ 

Fig.  18.  The  three  laminas  used  for  the  composite  laminate. 


Fig.  19.  Filtered  ^-acceleration  for  a  node  on  the  payload  using  a  0.2  in 
steel  plate  and  a  plate  with  0.2  in  thickness  and  30%  fiber  volume  fraction 
respectively. 


remainder  of  this  section  through  a  parametric  study.  In 
this  study,  the  fiber  volume  fraction  is  varied  from  30%  to 
70%  in  steps  of  10%  while  the  plate  overall  thickness  is 
varied  from  0.2  to  0.3  in.  in  steps  of  0.05  in.  Three  laminae 
of  equal  thickness  are  used  in  each  case.  Table  6  shows  the 
variation  of  the  material  properties  as  a  function  of  the 
fiber  volume  fraction.  The  table  shows  that  while  density 
and  Poisson’s  ratio  in  the  x-z  direction  experience  little 
variation,  the  values  of  Eyy,  and  are  more  than 
doubled  as  the  volume  fraction  increases. 

Finite  element  results  for  the  15  case  studies  are 
presented.  Table  7  shows  the  maximum  and  minimum 
filtered  axial  accelerations  outside  the  gun  barrel  for  the 
node  on  the  payload.  The  accelerations  on  the  payload  are 
reduced  significantly  using  a  composite  plate  when 
compared  to  the  original  steel  plate.  Fig.  19.  The 
acceleration  values  vary  over  the  considered  ranges  of 
composite  thickness  and  fiber  volume  fraction.  No  clear 
pattern  is  observed  to  recommend  specific  composite 
characteristics  over  others.  The  best  results  are  however 
obtained  from  the  following  cases: 

-  0.2  in.  thickness  and  70%  fiber  volume  fraction  (1753/ 

-1600g)  and, 

-  0.3  in.  thickness  and  60%  fiber  volume  fraction  (1651/ 
-1765  g). 

The  two  above  cases  correspond  to  cases  that  may  be 
difficult  to  make  and  maintain  control  of  quality.  If  such 
constraints  exist,  it  may  be  recommended  to  look  at 
alternatives  in  Table  7. 

Table  8  lists  the  frequency  of  the  filtered  axial  accelera¬ 
tion  for  the  15  case  studies.  Results  show  that  frequency 
increases  as  either  plate  thickness  or  fiber  volume  fraction 
increases,  which  is  consistent  with  the  change  of  the  plate 
mechanical  characteristics,  as  shown  in  Table  6,  and  the 
increase  of  the  bending  stiffness  of  the  plate,  respectively. 
Table  8  shows  over  60%  of  variation  in  frequency  range 
between  the  two  extreme  cases.  The  lowest  frequencies  are 
obtained  from  the  following  cases: 


0.2  in.  thickness 

and 

30% 

fiber 

volume 

fraction 

(1267  Hz)  and, 

0.2  in.  thickness 

and 

40% 

fiber 

volume 

fraction 

(1400  Hz). 


Table  6 

Variation  of  the  plate  material  properties  with  fiber  volume  fraction 


Fiber  volume  fraction  (%) 

Exx  =  E,,  (psi) 

Eyy  (psi) 

p  (Ibf/in^) 

Gxz  (psi) 

Gyy  (psi) 

Vxz 

30 

3.8E  +  06 

7.4E  +  05 

0.0502 

1.4E  +  06 

1.9E  +  05 

0.023 

0.33 

40 

4.9E  +  06 

8.6E  +  05 

0.0502 

1.8E  +  06 

1.9E  +  05 

0.019 

0.33 

50 

6.0E  +  06 

l.OE  +  06 

0.0540 

2.3E  +  06 

1.9E  +  05 

0.017 

0.32 

60 

7.2E  +  06 

1.3E  +  06 

0.0540 

2.7E  +  06 

1.9E  +  05 

0.017 

0.32 

70 

8.5E  +  06 

1.7E  +  06 

0.0579 

3.2E  +  06 

1.9E  +  05 

0.018 

0.32 
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Table  7 

Filtered  maximum  and  minimum  axial  accelerations  (G)  outside  the  gun 
barrel 


Plate  thickness  (in.) 

Fiber  volume  fraction  (%) 

30 

40 

50 

60 

70 

0.2 

1869 

-1910 

2109 

-2229 

2023 

-2061 

2350 

-2307 

1753 

-1600 

0.25 

1837 

-1719 

19449 

-19429 

1883 

-1821 

2166 

-2265 

1964 

-1575 

0.3 

2123 

-2076 

19719 

-1989 

2190 

-2153 

1651 

-1765 

2431 

-2340 

Table  8 

Frequency  of  filtered  accelerations  (Hz)  outside  the  gun  barrel 
Plate  thickness  (in.)  Fiber  volume  fraction  (%) 


30 

40 

50 

60 

70 

0.2 

1267 

1400 

1433 

1500 

1533 

0.25 

1600 

1667 

1733 

1767 

1800 

0.3 

1833 

1900 

1933 

2000 

2067 

Fig.  20.  Model  used  for  modal  analysis  of  the  payload  and  the  plate 
(0.2  in  thickness  and  30%  fiber  volume  fraction)  only. 


The  above  results  can  lead  to  the  recommendation  that 
relatively  thin  composite  plates  with  lower  natural  fre¬ 
quency  may  be  a  more  attractive  alternative  as  long  as 
stresses  are  within  design  limits.  To  verify  that  the  above 
results  depend  on  the  flexibility  of  the  gun  as  well  as  the 
mass  and  geometry  of  the  projectile,  a  modal  analysis  is 
conducted  for  the  payload  and  the  plate  (0.2  in.  thickness 
and  30%  fiber  volume  fraction)  only.  The  outer  edges  of 
the  plate  are  completely  fixed  as  shown  in  Fig.  20.  The  first 
ten  natural  frequencies  of  this  system  are  listed  in  Table  9. 
Comparison  of  the  results  show  that  modal  results  of 
simulating  the  gun  and  projectile  are  significantly  below 


Table  9 

Natural  frequencies  of  the  payload  and  the  plate  (0.2  in.  thickness  and 
30%  fiber  volume  fraction)  only 


Mode  number 

Frequency  (Hz) 

1 

2049 

2 

2368 

3 

2568 

4 

8181 

5 

8916 

6 

10,439 

7 

33,648 

8 

33,808 

9 

36,377 

10 

36,826 

Fig.  21 .  Mode  shape  for  the  first  mode  of  the  payload  and  the  plate  (0.2  in 
thickness  and  30%  fiber  volume  fraction)  only. 


those  of  the  plate  and  the  payload  only.  Fig.  21  shows  the 
mode  shape  of  the  first  frequency  of  the  plate  and  payload 
(Table  10). 

5.  Conclusions 

Components  within  a  projectile  are  subject  to  severe 
shock  loadings  during  and  immediately  after  a  gun  launch. 
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Table  10 

Equivalent  SI  units 


US  units 

SI  units  equivalent 

Length 

1  inch 

0.0254  m 

Weight 

1  Pound 

0.4535  kg 

Temperature 

IF 

(5/9)*(lF-32)°C 

Acceleration 

1  in/s^ 

0.0254  m/s^ 

Specific  weight 

1  Ibr/in^* 

273386.3  N/m^ 

Young’s  modulus 

1  psi 

6,894.75  N/m^ 

Density 

1  Ib/in^ 

27,679.9  kg/m^ 

Extensive  experimental  tests  of  such  components  can  be 
prohibitively  expensive.  This  paper  presents  an  approach 
to  model  the  gun-launch  dynamic  interaction  between  a 
projectile  body,  its  internal  components,  and  the  gun 
barrel.  The  proposed  approach  includes  the  effects  of  the 
flexibility  of  the  gun  barrel,  friction  between  the  projectile 
and  the  gun,  and  flexibility  of  a  mounting  plate  for  an 
electronic  payload.  While  FEM  displacement  and  velocity 
results  for  nodes  on  projectile  were  satisfactory,  there  is  a 
need  for  filtering  the  FEM  acceleration  data  as  the 
acceleration  results  contained  higher  level  of  noise  than 
what  is  observed  in  corresponding  test  data.  A  filtering 
approach  based  on  modal  analysis  of  the  projectile  is  used 
in  this  research. 

A  parametric  study  is  also  conducted  to  explore  the 
possibility  of  using  a  composite  plate  to  support  an 
electronic  payload.  In  this  study,  fiber  volume  fraction 
and  thickness  are  varied.  All  cases  in  the  parametric  study 
experience  marked  reduction  in  magnitude  and  frequency 
of  payload  accelerations  when  compared  to  the  results 
found  when  using  a  steel  plate  of  the  same  size.  Simulation 
results  also  show  that  the  composite  plates  with  lower  fiber 
volume  fraction  and  thickness  provide  the  lowest  frequency 
of  the  acceleration  response  after  the  gun  muzzle  exit.  This 
observation  can  be  useful  in  reducing  the  possibility  of  low- 
cycle  fatigue  failure  of  electronic  components. 

The  current  model  is  computationally  intensive.  Meth¬ 
ods  for  reduced  computational  time,  without  sacrificing 
accuracy,  are  currently  under  investigation.  Experimental 
verification  of  shock  transmission  through  composite 
plates  is  needed  for  further  tuning  of  the  FEMs.  Future 
studies  may  include  effects  of  gun  recoil  on  the  model  and 
fiber  orientation  angles  of  the  laminae. 
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Appendix  A.  Elastic  properties  of  composite  laminates 


Elastic  properties  of  the  composite  lamina  are  estimated 
using  Rule  of  Mixture’s  equations  following  [17].  Young’s 
Modulus  of  the  composite  lamina  along  the  fiber  direction 
is  calculated  using  the  following  equation: 

En  =  EfVf  ^  (A.l) 

Transverse  modulus  and  shear  modulus  of  the  composite 
lamina  are  calculated  as  follows: 


„  EfEra 

EfV^  +  E^Vf’ 

(A.2) 

^  GfGm 

GfV^  +  G^Vf- 

(A.3) 

Major  and  minor  Poisson’s  ratios  of  the  composite  lamina 
are  calculated  using  the  following  formulas: 

^12  =  Uf  “h  U m?  (^-4) 


f^2i  =  ^  Vi2.  (A.5) 

-fcll 

Fig.  18  shows  the  stacking  sequence  of  the  laminas  in  the 
laminate.  The  material  properties  of  the  laminate  are  found 
by  calculating  extensional  stiffness  matrix  [A]  of  the 
laminate. 


[A]  = 


11 

^12 

^16 

12 

^22 

^26 

16 

^26 

^66 

(A.6) 


[A]  is  calculated  using  the  following  formula: 

[A]  =  [2]+6oo(^o  -  hi)  +  [Q]Qo(hi  -  h2)  +  [2]_6o°(^2  -  h). 

(A.7) 

where  is  the  stiffness  matrix  for  the  individual  lamina, 
which  is  defined  as 


[Q\e  = 


'll 

2i2 

Q\6 

'12 

2i2 

Q26 

16 

Qi6 

066 

(A.8) 


Elements  of  [Q\q  can  be  calculated  using  the  following 
formulas  [17]: 

Qn  =Qn  cos'^  0  +  2(gi2  +  2g66)sin^  0cos^  0  +  Q22  sin'^  0, 

(A.9) 

612  =  Qi2(sir2  0  +  cos^^  0)  +  (fill  +  222  -  4266)sin^  0  cos^  0, 

(A.  10) 

0.22  =  On  0  +  2(2i2  +  2266)sin^  0cos^  0  +  222  0, 

(AM) 
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Qi6  =  (2ii  -  Qn  -  2266)  sin  0cos^  9 

+  (2i2  -  222  +  2266)sin^  0  cos  0,  (A.12) 

226  =  (2ll  ~  2i2  ~  2266)sin^  0  cos  6 

+  (2i2  -  222  +  2266)  sin  0COS^  0,  (A.13) 

266  =  (2ll  +  2i2  -  22i2  -  2266)sin^  0cos^  e 

+  266(sin'^  0  +  cos"^  0),  (A.  14) 

where  0  is  the  fiber  orientation  angle  of  the  lamina 
En 


Qn  = 

Qii  = 


1  —  ^9121^21  ’ 

Ell 

1  -  DiiDii  ’ 

^iiEii 


1  —  ^9121^21 


2i2  —  621  — 

Qee  =  ^12- 

The  properties  of  the  laminate  are: 
-  In-plane  Young’s  modulus: 


Er.x  =  E,,=- 


^11 


(ho  +  hi)  All 

-  In-plane  Poisson’s  ratio: 

^12 

^11 

-  In-plane  shear  modulus: 

^  _  ^11  -^12 

2{ho  hi) 

-  Transverse  modulus: 

Eyy  =  Ell. 

-  Poisson’s  ratio  in  the  transverse  direction: 

^xj;  =  1^21. 

-  Shear  modulus  in  the  transverse  direction: 


^xy  — 


Density  of  the  Composite  Laminate: 


(A.15) 


(A.  16) 


(A.  17) 


(A.  18) 


(A.  19) 


(A.20) 


(A.21) 
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